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Abstract 

Qh We point out that existing numerical data on the correlation length 

and magnetic susceptibility suggest that the two dimensional 0(3) model 
with standard action has critical exponent tj = 1/4, which is inconsis- 
tent with asymptotic freedom. This value of tj is also different from 
the one of the Wess-Zumino-Novikov-Witten model that is supposed to 
correspond to the 0(3) model at 9 = it. 



Ever since the seminal paper of Kosterlitz and Thouless [1] it has been 
accepted that in two dimensions (2D) there is a fundamental difference be- 
tween the Abelian 0(2) and the non-Abelian 0(3) models: while the hrst one 
undergoes a transition (KT) from a phase with exponential decay at high tem- 
perature T(T = 1/(3) to a massless one at low temperature, the latter is in the 
high temperature phase for any fi < oo. While the hrst conjecture was proven 
rigorously by Frohlich and Spencer [2], the second one remains unproven. 

We have advanced several arguments [3, 4, 5] against this accepted belief 
and in our opinion all 0(N) models undergo a KT-like transition at some finite 
(3. Since however our arguments do not constitute a rigorous demonstration, it 
may be useful to get a glimpse at what the truth may be via existing numerical 
data. Our observation pertains to the value of the critical exponent t], defined 
as follows: let £ be the correlation length and x the magnetic susceptibility. 
Then r/ can be defined by the following asymptotic statement: 
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i.e. In x — (2 — Tf) In £ should go to a constant as £ — > oo. 
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Figure f: 0(2) model, test of i] = 1/4 
For the 0(2) model certain non-rigorous renormalization group arguments 
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advanced by Kosterlitz [6] predict that r/ = 1/4. In Tab.l we collect the data 
from [7] and [8]. We could plot ln% vs. ln£ and would see a nice straight line 
with a slope close to 1.75. In fact a fit of In % to a linear function of ln£ yields 
a slope of 2 — r/ = 1.72 with a good fit quality, if the errors on both £ and x 
are taken into account. It is, however, much more significant to look at the 
deviations from the predicted slope by plotting In % — 1.75 In ^ vs. ln£, together 
with the least square fit to a constant. This is done in Fig.l; the data of [7] are 
represented by crosses, the ones of [8] by squares. It can be seen that the data 
still show a tendency to decrease, corresponding to the fact that the aymptotic 
region has not yet been reached and the fits published in the literature always 
produced an r/ > 1/4, as did ours. But from Fig.l one can also see a tendency 
of the data to flatten with increasing £, in accordance with the theoretical 
prediction of r/ = 1/4. The error bars shown in Fig.l correspond to adding the 
errors on ln% to 1.75 times the errors on ln% (probably an overestimation as 
the errors on x an d £ are correlated). 

According to the accepted scenario, the behavior of the 0(3) model should 
be vastly different; indeed perturbative renormalization group arguments pre- 
dict the following asymptotic behavior ("asymptotic scaling"): 



From this it follows that for £ — > oo one should observe In x — 2 In £ + 2 In (3 
tending to a constant. To test this prediction, we use the published data of 
[9] as well as the unpublished ones due to [10], which are collected in Tab. 2. 
In Fig. 2 we plot In x — 2 In £ + 2 In /3, which should be constant according to 
the asymptotic scaling prediction. In Fig. 3 we plot instead ln% — 1.75 In £ vs. 
ln£, as we did for 0(2). Fig. 3 also contains a least square fit to a constant. 
The data from [9] are represented by squares, the ones of [10] by crosses. The 
error bars are obtained, as in the case of 0(3), by adding the erros on ln% and 
1.75 times the errors on ln£. Note that in both figures we are using the same 
scale as in Fig.l. Fitting In % to a linear function of ln£, as we did for 0(2), 
gives a slope of 2 — r/ = 1.74, again with a good fit quality (% 2 1 per degree 
of freedom). 

The following facts are visible from Figs 3 and 2: 

(1) The data are not consistent with asymptotic scaling inasmuch the points 
shown in Fig. 2 are not consistent with a constant but show a strong decrease. 

(2) The data are consistent with a critical behavior with the same r/ = 1/4 
that was predicted by the KT theory for the 0(2) model. 

So we have to conclude that in fact the numerical data disagree with the 
predictions of asymptotic freedom and instead suggest that both 0(2) and 
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Figure 2: 0(3) model, test of asymptotic scaling 



(9(3) have the same value of r/. The data also indicate that the 0(3) model 
with standard action J2 s (i) ' s (j) 1S n °t in the same universality class as the 
6 = 7r 0(3) model, which is supposed to have the Wess-Zumino-Novikov- 
Witten model as its scaling limit [11]; conformal held theory arguments predict 
J]=l, a predicition which a recent paper [12] claims to have verified numerically. 

A note of caution: in principle r/ can also be determined by studying the 
finite size scaling of the susceptibility x(L)- Ideally this study should be per- 
formed at (3 crt ; in practice though, it suffices to place oneself in a regime where 
£ >> L. One could ask then what value of r/ would come out of such a de- 
termination. The answer is that one must be careful and use a large enough 
L. Indeed, in this type of problems, besides the correlation length there 
is a second important length £px, the distance over which the system is well 
ordered. From perturbation theory (PT) one knows that for the 0(3) model 
this second length is 0(exp(7r/3)). Consequently, if one uses periodic boundary 
conditions (b.c.) and an L and (3 such that ^pt >> L, then the system will 
behave according to PT and one will find the perturbative value of r/ which 
is I/tt/3 + 0(/3~ 2 ). As we explained in a recent paper [14], a clear indication 
that such a determination cannot be trusted to reflect the true thermodynamic 
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Figure 3: 0(3) model, test of i] — 1/4 



value of ?7 is the dependence of r/ upon the b.c. employed. Indeed, as we dis- 
cussed in ref. [14], in the limit L — > oo, the value of r/ must be the same if we 
use periodic b.c. and freeze two spins situated a distance L/2 apart in any 
arbitrary positions; this then gives a criterion for deciding whether one has 
used a sufficiently large L and whether the determination of r/ is trustworthy. 

Finally we should like to make a remark about the large iV-limit of the 
O(N) model, which is often cited as evidence for the standard picture of (3 crt = 
oo for all N > 2. As we pointed out in various places (see in particular [15]), 
this reasoning is faulty. It is true that at fixed (3 = (3jN both £ and x converge 
as N — y oo to the values of the spherical model, hence in that limit r/ = up 
to log corrections. But, as we pointed out in [15], the approach to the limit 
is nonuniform, i.e. it becomes slower and slower as (3 — > oo or £ — > oo. 
With increasing N, In % as a function of ln^ will approach the spherical limit 
(iV = oo) form, but at any fixed N, for large enough f3, we expect that it will 
fall away from that form and that ln% — 1.75 In £ will aproach a constant for 
(3 — > oo. 

We are indebted to S.Caracciolo, R.G.Edwards, A.Pelissetto and A.D.Sokal 
for allowing us to use their data prior to publication; these data form the basis 



4 



of the papers [13] and will apear in full in subsequent papers by [16]. 
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Tab. la: x an d £ f or the 0(2) model on thermodynamic lattices (L/£ > 7). 
Dat taken from [7]. 



p 


f 

S 


y 


.91 


9.36(5) 


134.17(45) 


.92 


10.69(8) 


162.68(61) 


.93 


12.03(6) 


201.29(53) 


.94 


13 50(9) 


249 90(80) 


.95 


15.61(10) 


319.6(1.1) 


.96 


18.08(13) 


414.1(1.7) 


.97 


21.66(13) 


554.9(2.9) 


.98 


26.37(19) 


764.6(3.0) 


.99 


31.78(21) 


1092.(4.) 


1.00 


40.20(20) 


1604.(4.) 


1.02 


69.27(59) 


4170.(19.) 



Tab. lb: x an( t £ f or the 0(2) model on thermodynamic lattices (L/£ > 7). 
Dat taken from [8]. 



p 




X 


.9090909 


9.32(2) 


132.41(55) 


.9345794 


12.61(2) 


220.84(86) 


.9478673 


15.23(5) 


302.19(89) 


.9615385 


18.70(20) 


433.7(2.2) 


.9708738 


21.80(20) 


564.1(2.7) 


.9756098 


23.90(20) 


665.0(4.0) 


.9803922 


26.20(20) 


779.0(5.0) 


.9900990 


31.70(30) 


1087.(10.) 


1.000000 


40.40(40) 


1631.(12.) 


1.010101 


51.30(90) 


2520.(24.) 


1.020408 


70.(1.) 


4258.(20.) 



Tab. 2a: x an( t £ f or the 0(3) model on thermodynamic lattices (L/£ > 7). 
Dat 3 31*6 taken from [9]. 



p 




X 


1.5 


11.05(1) 


176.4(2) 


1.6 


19.00(2) 


448.4(7) 


1.7 


34.44(6) 


1263.7(3.3) 


1.75 


47.2(2) 


2197.(15.) 


1.8 


64.5(5) 


3823.(21.) 


1.85 


88.7(5) 


6732.(25.) 


1.9 


122.7(1.1) 


11867.(62.) 
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Tab. 2a: x an d £ f or the 0(3) model on thermodynamic lattices (L/£ > 7). 
Data are taken from [10]. 



1 




/V 


1.5 


11.13(16) 


178.0(1.1) 


1.55 


14.69(15) 


281.5(1.7) 


1.6 


18.66(35) 


447.6(2.6) 


1.65 


25.21(13) 


742.8(2.0) 


1.675 


29.77(13) 


971.9(2.7) 


1.7 


34.66(13) 


1269.4(3.6) 


1.725 


40.21(25) 


1661.8(5.8) 


1.75 


47.20(19) 


2193.9(6.0) 


1.775 


55.23(25) 


2886.4(10.2) 


1.8 


65.17(57) 


3839.7(23.8) 
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